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Technologies are being developed to manipulate electromagnetic waves using artificially structured 
materials such as photonic crystals and metamaterials, with the goal of creating primary optical devices. For 
example, artificial metallic periodic structures show potential for the construction of devices operating in 
the terahertz frequency regime. Here we demonstrate the fabrication of photo-designed terahertz devices 
that enable the real-time, wide-range frequency modulation of terahertz electromagnetic waves. These 
devices are comprised of a photo-induced, planar periodic-conductive structure formed by the irradiation of 
a silicon surface using a spatially modulated, femtosecond optical pulsed laser. We also show that the 
modulation frequency can be tuned by the structural periodicity, but is hardly affected by the excitation 
power of the optical pump pulse. We expect that our findings will pave the way for the construction of 
all- optical compact operating devices, such as optical integrated circuits, thereby eliminating the need for 
materials fabrication processes. 

Terahertz radiation (1 THz = 10 12 Hz), which is commonly defined as the region of the electromagnetic 
spectrum between 0.3 THz and 10 THz, is extremely attractive for a variety of applications in science and 
technology such as sensing, imaging, and telecommunication 1 . Because exploration of this spectral region 
located between radio and optical frequencies is only at the initial stages, many further applications in diverse 
fields are also expected. One envisaged technology involves the development of high-speed and high-frequency 
operating devices in the field of ultrafast optical communication, which has the potential to replace conventional 
electric circuits. In order to realize applications such as optical switching devices, control of the frequency 
characteristics over a wide frequency range is important. 

Previous studies have verified that THz pulses can be guided using metal wires 23 , and that the frequency 
modulation can be tuned by artificial periodic metal structures 4 ' 5 when the geometry and dimensions are appro- 
priately scaled 6 ' 7 . Furthermore, the real-time optical control of THz pulses has been achieved by the application of 
a bias voltage in so-called active terahertz metamaterials 8 , by the injection of photo-induced charge carriers into 
metamaterial arrays 910 , and by the use of reconfigurable metamaterials that respond to external stimuli 11 . The 
dynamic control of metamaterials has great potential for use in applications such as optical modulators with 
temporally and geometrically varying patterns, as well as in the active control of optical filters 1213 . However, the 
modulation frequency has thus far been limited to a small region around the resonance frequency determined by 
the fixed structural geometry, allowing only narrow frequency shifts. In order to overcome this limitation, we 
propose and demonstrate the operation of a photo-designed THz device in which the period, shape and dimen- 
sion of the photo-induced planar structural pattern can be varied dynamically. 

Results 

The planar photo-induced structures used in this study were one-dimensional periodical rectangular corruga- 
tions (namely, line grating (LG) patterns), which were created by focusing a spatially phase-modulated femto- 
second laser pulse on a high-resistivity Si (R > 1000 Q»cm) surface. The fabrication system and the geometry of 
the LG structure are shown in Fig. la. The optical pump pulse wavelength of 800 nm allows the excitation of 
carriers to the conduction band, the spatial distribution of which on the silicon (Si) surface can be controlled using 
a liquid- crystal spatial light modulator (SLM), as shown in Fig. lb (see Methods). This enables the creation of 
arbitrary beam patterns with periodicities that can be varied freely. The lines comprising the LG pattern were well 
defined, as shown in the optical microscope image in Fig. lc, although small spatial fluctuations of the photo- 
induced carrier density might be inherent 14 . The period of the LG structure was varied from a = 100 um to 
300 um; the width of each line in Fig. lc corresponds to one-half of this period. The optical pump power (P pump ) 
was varied from 4 mW to 100 mW, which gave a corresponding carrier density at the Si surface (iV pump ) of 
between 1.1 X 10 16 cm" 3 and 2.7 X 10 17 cm" 3 , and a plasma frequency (f p ) of between 0.7 THz and 3.4 THz. 
These values were calculated using the relationship f p = (iV pump e 2 /s 0 e Si u) 1/2 , where l/|i = l/m e + l/m h is a 
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Figure 1 | Fabrication and experimental operation of photo-designed THz devices, (a) Fabrication of a planar photo-induced LG pattern at the total 
reflection surface of a Si prism with refractive index 3.4 in the THz frequency regime, (b) Measurement of surface electromagnetic mode on photo- 
designed THz devices (see Methods), (c) Optical microscope image of photo-induced LG pattern taken with an infrared sensor viewer for better contrast. 
Fluctuations in brightness within the lines are due to scattered light at the rough surface of the viewer. The periodicity of the pattern is parallel to the THz 
wave vector fc THz . 



function of the reduced masses of the electron and the hole, 
m e = 0.27m 0 and m h = 0.37m 0 are the effective masses of electrons 
and holes in Si 15 , s 0 is the permittivity of free space, s Si = 11.6 is the 
permittivity of Si in the THz frequency regime, e is the electron 
charge, and m 0 is the electron rest mass. 

We employed an internal-reflection configuration for our mea- 
surement (see Methods), checking that no undesired THz emission 
from optical excitation inside the Si occurred 16 . The optical-pump 
and THz-probe technique was used to characterize the performance 
of the photo-designed THz devices, as shown in Fig. lb. We note that 
diffraction between the pump and probe pulses, as observed for 
holograms and population gratings 1417 " 19 , should be prevented here 
because the frequency of the optical pulses was much higher 
than that of the probe pulses. The timing of the optical excitation 
was At = 50 ps before the arrival of the THz probe pulse at the Si 
surface; this time delay was large compared with the rise time of the 
photocarriers, and sufficiently shorter than the photocarrier lifetime 
in bulk Si crystals 20 ' 21 . Furthermore, diffusion effects were expected to 
be negligible, because the diffusion length of the photocarriers is 
0.5 |im at 50 ps after excitation 19 . Therefore, the quasi-steady-state 
properties of the photo-designed THz devices could be investigated 
as a function of the carrier density, which was determined by the 
excitation power. Reflectance spectra were obtained using THz time- 
domain spectroscopy, which is described elsewhere 1,22 . 

The frequency- dependent reflected intensity is plotted as a func- 
tion of the LG period in Fig. 2. Here, P pU mp was set to a constant value 
of 40 mW, for which the corresponding carrier density at the Si 
surface is N pump ~ 1.1 X 10 17 cm" 3 and the corresponding plasma 
frequency isf p — 2.1 THz. No characteristic features, except for a 
broad minimum around 2 THz, were observed in reflection spec- 
trum of a reference photo-induced film (in which the whole area 
occupied by the LG pattern was uniformly irradiated at the same 
carrier density of AT pump ~ 1.1 X 10 17 cm" 3 ). By contrast, a sharp 
minimum at lower frequencies was observed for all the LG structures; 
no broad minimum was evident. No significant changes in the spec- 
tra were observed upon varying At between 5 ps and 300 ps. This 
was consistent with previously reported ambipolar-diffusivity mea- 
surements made using the population-grating method, from which a 
carrier diffusion length of 1.3 |im was estimated at 300 ps (ref. 19). 



As shown in Fig. 2, the frequency of the characteristic sharp 
minimum could be widely modulated from 0.4 THz to 1.2 THz by 
deliberately decreasing the period of the LG pattern. 

The dependence of the modulated frequency on the periodicity is 
shown in Fig. 3. The frequency of the sharp minimum increased as 
the periodicity decreased, in good agreement with theoretical con- 
siderations, where the surface electromagnetic mode is assumed to be 
on the periodic structure 23 " 26 . The theoretical curve in Fig. 3 was 
obtained using the approximation from the modal expansion 
method 25 , where / sw ~ c/[a(8 Si 1/2 sin 0 +1)], under the reasonable 
assumption of a thin, conductive and planar LG structure; c is the 
velocity of light in vacuum and 0 is the incident angle of 30°. Our 
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Figure 2 | Modulation performance of photo -designed THz devices as a 
function of the period of the photo-induced LG. The solid lines represent 
a series of reflectance spectra from photo-induced LG structures with 
various periods. The dotted line represents the spectrum for the reference 
photo-induced film at the same carrier density. The reflection spectrum 
with no optical pump pulse irradiation is employed in order to normalize 
the data. 
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Figure 3 | Dependence of spectral minimum frequency on structural 
periodicity in photo-designed THz devices. Schematic representations of 
the one-dimensional LG (blue circular data points) and the two- 
dimensional square array patterns (red square data points) are shown in 
the inset. White and black areas indicate irradiated and non-irradiated 
parts, respectively; the wave vector of the THz pulses fc THz is in the 
horizontal direction. The theoretical curve obtained under the assumption 
of the surface electromagnetic mode is shown as a dashed line. 

results thus imply that the photo-induced structure is similar to a 
planar conductive structure. We also investigated the dependence of 
the spectral minimum on the periodicity of a two-dimensional 
square array (each square in the illustration shown in the inset of 
Fig. 3 occupies half a period). The corresponding reflectance spectra 
are shown in Supplementary Fig. SI. The same trend of the char- 
acteristic spectral minimum with respect to periodicity was observed 
for the square array, suggesting that the frequency of the minimum is 
primarily dependent on the periodicity of the structure and not on 
the pattern, because the structure is much thinner than the THz 
wavelength due to the skin depth of Si (12.2 |im) at a wavelength 
of 800 nm 27 . 

Because the photo-induced carriers forming the LG structure 
occupy a smaller surface area interacting with the THz probe pulses 
than the reference photo -induced film, one would expect that the 
attenuation due to conductivity losses is lower than in the reference 
film 2 . Therefore, a sharp signal from the surface electromagnetic 
mode should become clearly observable. Nevertheless, the conduc- 
tivity losses are still inherently associated with the photo-induced 
structure 13 over the entire THz probe frequency range. It is inform- 
ative to consider the dependence of reflectance on the photon - 
induced carrier density. 

Figure 4 shows the transient evolution of the reflectance spectra at 
different values of pump photon power (P pump ) in the case of a 
photo-induced LG with a = 170 |im (the calculated frequency of 
the minimum is 0.65 THz). As the excitation pump power increased, 
the reflectance dip (the spectral minimum in the region of 0.7 THz) 
initially became deeper and sharper, with an accompanying decrease 
in reflectance over the entire frequency range. The spectral minimum 
was blue shifted from 0.54 THz to 0.83 THz as P pump increased, and 
the minimum value of reflectance was found at P pump — 50 mW 
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1.3 X 10 17 cm" 3 and/ p « 2.4 THz). With further increase 



in ^pump* the reflectance dip decreased and eventually became incon- 
sequential. We note that the spectral minimum appears only at 
P pump > 10 mW, a value corresponding to/ p ~ 1.1 THz, which is 
larger than the frequency of the spectral minimum (f p >/ S w)- This 
implies that the spectral minimum does not arise from modulation of 
the refractive index as occurs in photonic crystals and electromag- 
netically induced gratings 17 " 19 . 
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Figure 4 | Modulation performance of a photo-designed THz device as a 
function of optical pump power. The period of the LG structure measured 
was fixed at a = 170 um. A minimum frequency of 0.65 THz was 
calculated under the assumption of the surface electromagnetic mode. The 
reflection spectrum with no optical pump pulse irradiation is employed in 
order to normalize the data. The inset shows the modulation efficiency of 
the reflectance at the minimum as a function of optical pump power. The 
reference reflectance was set at a frequency of 1.5 THz, indicated by the 
black arrow. 

Discussion 

In order to interpret these observations quantitatively, we evaluated 
the modulation efficiency of the reflectance given by -log 10 |#e X (/min)/ 
^ex(/ref)| 2 > where R ex (f) is the reflectance, / min is the frequency at the 
spectral minimum, and/ ref is the reference frequency; here, we used 
1.5 THz (far enough from the minimum) as a typical f Te{ . The modu- 
lation efficiency reaches a maximum at P pum p ~ 50 mW, as shown in 
the inset of Fig. 4, corresponding to the power at which the sharpest 
minimum was observed. The frequency here is slightly higher than 
the calculated value of 0.65 THz. This discrepancy cannot be due to 
the finite thickness of the structure because the frequency of the 
minimum should be constant with respect to thickness 28 . In the 
internal reflection configuration used in our study, residual carriers 
excited by the previous pump pulse (1 ms before the current mea- 
surement) may disturb the reflectance spectrum under high-power 
excitation (Supplementary Figure S2). The blue shift of the minimum 
would then be caused by a reduction of the effective permittivity of Si 
due to the existence of diffused residual carriers. Because the density 
of diffused residual carriers depends on the photo -irradiated surface 
area, the excitation intensity should be optimized for each particular 
photo-induced structure in order to achieve high modulation effi- 
ciencies for practical applications. We note that the frequency of the 
minimum can be significantly modulated by the period of the struc- 
ture, but the excitation pump power has a much smaller influence. 
This suggests that the sharp minimum cannot be attributed to the 
waveguide mode that characteristically shifts on varying the thick- 
ness of the structures 28 . 

For the realization of higher performance photo-designed THz 
switching devices, the quality factor of the surface electromagnetic 
mode (the width of the minimum) should be improved. One con- 
ceivable approach is to reduce the surface area of the photo - 
irradiated lines without changing the period. In addition, further 
optimization of the spatial distribution of the photo-induced carrier 
density in the photo-irradiated part of the structure should also be 
effective, for example, by the utilization of computer-generated holo- 
grams calculated using an iterative Fourier transform algorithm 14 . It 
is important to note that photo- designed THz devices would also be 
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usable in a transmission configuration with normal THz incidence. 
In this configuration, modifications to control the thickness of the 
structure that interacts with the THz probe pulses would be necessary 
in order to increase the modulation efficiency, for example, two- 
photon excitation. This solution would be applicable to our experi- 
mental configuration because the densities of the photo -induced 
carriers and the corresponding plasma frequencies exponentially 
decrease with depth inside the Si prism in the direction perpendicular 
to the surface 22 . Moreover, the maximum modulation speed in 
switching devices is limited by the lifetime of the photo -induced 
carriers in the Si or the response time of the SLM. The bulk recom- 
bination lifetime of the carriers is of the order of 10 jas 20,21 at a low 
optical pump power, which could improve the operating speed of 
devices such as present-day electrical switching devices 7 ' 29 . However, 
ultrashort lasers such as femtosecond pulse lasers are preferable and 
indispensable for practical use in photo -designed THz devices. Rapid 
and stable modulation performance using time and space variations 
of the photo -induced structures can be realized by brief optical 
excitation in order to prevent distribution changes of the photo - 
carrier 19 ' 22 during the modulation, while still generating sufficient 
photocarrier density. The switching speed could be improved by 
choosing a suitable semiconductor, for example a direct semi- 
conductor such as GaAs, or by using a thin Si film 30 , which would 
also reduce the conductivity losses. In the case of semiconductors 
with shorter bulk recombination lifetimes such as low-temperature- 
grown GaAs, it is also necessary to consider the limitation that arises 
from the greater surface recombination velocity, which causes a 
much faster decrease of the surface carrier density. The use of a 
ferroelectric liquid- crystal SLM might also allow the switching speed 
to be improved. Furthermore, a semiconductor with a small self- 
diffusion coefficient would be an appropriate choice. This solution 
could enable the utilization of an electrically modulated continuous - 
wave laser diode as the optical pump source in photo -designed 
THz devices. 

This work represents the first step towards realizing a novel tech- 
nology for the wide- range frequency modulation of the THz regime, 
as well as towards establishing techniques for the real-time operation 
of the photo -induced structures that make frequency modulation 
possible. That is, we have established a route towards the real-time 
creation, modulation, and removal of arbitrary photo -induced struc- 
tures while eliminating the need for fabrication processes or electrical 
connections. The experimental demonstration of our technique has 
other significant implications, such as the ability to control optical 
responses by dynamical breaking of the symmetry properties of the 
crystals. This work also has potential to open new technology to 
realize hybrid switching devices that utilize both conventional sur- 
face plasmons with flat surfaces and mimicking (or designer) surface 
plasmons 24 ' 26 with structured surfaces. 

Methods 

The surface electromagnetic mode of our photo -designed THz devices was investi- 
gated using the internal- reflection configuration (Fig. lb) withp-polarized THz probe 
pulses (incidence angle 30°). This setup uses a non-coaxial configuration between the 
pump and probe pulses, which is convenient for both practical use and better 
modulation efficiency. Propagation of the surface electromagnetic mode gives rise to 
a large modulation of the THz reflection signals. The modulation efficiency can be 
quantitatively investigated as a function of frequency, because the intensity of the THz 
incident wave can easily be obtained by cutting the optical pump pulse (the attenuated 
total reflection configuration), which is not possible using the normal reflection 
configuration. The beam waist of the THz pulses at the Si surface was 3.8 mm (smaller 
than the photo -irradiated area), and the Rayleigh range was 76.4 mm, which indicates 
that the THz pulse can be regarded as a plane wave (these values were calculated for a 
frequency of 1 THz). The THz probe wave vector /cthz (transverse magnetic or TM 
mode) was perpendicular to the LG structure (Fig. lc). The optical pump pulses were 
generated by a Ti:sapphire regenerative amplifier (Spitfire, Spectra Physics Inc.) with 
a centre of wavelength of 800 nm, a power of 800 uj/pulse, a repetition rate of 1 KHz 
and a pulse duration of 130 fs. The same system was used to generate the p-polarized 
THz probe pulses. The beam pattern focused at the Si surface to create both the photo- 
induced LG and square array structures was modulated using a phase-only liquid- 
crystal SLM (PLUTO, HOLOEYE Inc.). The computer-generated holograms of the 



LG (or square array) patterns were obtained using an inverse 2D Fourier transform 
algorithm, addressed to the SLM, and then stored on it. The phase of the laser pulse 
from a regenerative amplifier is spatially modulated by the reflection at the SLM. The 
SLM is separated from the lens shown in Fig. lb, a so-called Fourier lens, by the focal 
length of the lens. The 2D Fourier transform of the phase-modulated pulse from the 
SLM is performed by the Fourier lens, yielding the LG (or square array) intensity 
patterns on the Si surface. A reference film, in which the whole area occupied by the 
LG pattern was uniformly irradiated, was also examined. In this case, the SLM was 
replaced by a mirror, and the lens shown in Fig. lb was removed. 
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